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1. INTRODUCTION 

In today’s world, it’s always been a concern for an engineer to get a product out at a very optimal 
cost by minimizing both the product operating cost and raw material input to the production unit. Economic 
Load Dispatch (ELD) deals with the same situation and it works on operating a coordinated power system 
such that the lowest operating cost generators are used to the greatest extent and the highest operating cost 
generator is used to the lowest extent. Economic load dispatch problem is a constrained problem, several 
principles and strategies are already being developed to solve these problems. ELD has become an important 
fundamental function in operation and control of the power system. 

The demand for electricity is increasing in a large factor in today’s life, which makes it highly 
crucial to run generators at very minimal cost. This is the main factor of an Economic dispatch problem. With 
the unexceptional production of carbon emissions in thermal power plant, its needed to optimize the emission 
together with the optimization of cost which acts as two vital parts of Economic dispatch problem. The 
economic dispatch solution provides the best minimum cost of fuel and emission. This indirectly makes 
lower cost for electricity and makes electrical utilities more competitive in the market. As the energy cannot 
be stored, it requires highly efficient estimation scenarios including transmission and distribution systems to 
make the same work effectively. 

Various technologies have been introduced to solve the optimization of Economic Load Dispatch 
problems. The selection of the optimization algorithm is the important part of the problem involving 
economic dispatch. The EDP is developed based on real-valued codification. In modern methodology only 
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the cost function is evaluated and a global minimum solution is computed, independently of the cost function. 
The use of digital computers for obtaining loading schedules were investigated and used today. 

Many deterministic optimization approaches were proposed to solve the ELD problem, including 
lambda iteration method [1], gradient method, linear programming [4], non-linear programming, dynamic 
programming [2] and quadratic programming [14]. But these methods require enormous efforts in terms of 
computation. Due to complexities of computing, therefore efficient algorithm to find optimal solution like 
genetic algorithm [16], [18], particle swarm optimization [5], evolutionary programming, artificial bee 
colony optimization [9], [10], and biogeography based optimization; bacterial foraging and also their variants 
came into implement. Bio-inspired meta-heuristic algorithms have recently shown the efficiency in dealing with 
many nonlinear optimizations constrained problems for finding the optimal solution. 

Recently a nature based optimization technique whale optimization Algorithm (WOA) is developed 
based on the flashing behavior of Whales. WOA developed and is used to solve constrained engineering 
problems. Until now many researches have been carried out to find the closest optimum result in determining 
the power generation of each generator using WOA and it was inferred that the WOA is more robust and 
efficient in determining the optimal load scheduling. 


2. PROBLEM FORMULATION 
The generating units are loaded economically such a way to reduce the operating cost. Considering 
the valve point effect the economic dispatch formulated the objective function as given below 


FC; = a; + bP; + ey (1) 
Where 


Gis Bis Ci are the fuel cost coefficients of generator’ 
F; is the power generated by unit *, MW 


F Ci is the fuel cost function of unit * 


The total fuel cost for the entire system of N generators can then be calculated as, 


N N 
FC, = ». FGP) = >. a; + b,P, + cP? 
=1 i=1 (2) 
The new objective function by considering valve point loading along with total fuel cost becomes, 
N 
FC; = > a;+ bP, +c + |d; sin( e; ( P:min — P ))| 
=1 (3) 


Where 


a;, by Cis d; and ©; are the fuel cost coefficients of generator i 
P% is the power generated by unit È, MW 
Ni is the number of generating units 
Pi minis the minimum generation limit of unit È, MW 
F C; is the fuel cost function of unit £ 
F Cy is the total fuel cost, $/hr. 
In order to minimize the pollutants, emission is consiered along with economic dispatch. The generator can 


be modelled as having a quadratic relation between the amount of pollutants released and the power 
generated. 
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The mathematical formulation for generator is given by, 
— 2 
E; = a, + BP, Hy: 
Where 


Zin Bis Yi are the emission coefficients of generator © 
P; is the power generated by unit *, MW 


E; is the fuel cost function of unit = 


The total emission for the entire system of N generators can then be calculated as, 


N N 
Er = > EŒ) = > a, + BP, +y;,P? 
1 1 


The new Emission function becomes, 


N 
Er= Y a; + BiP, +y:P? + eexp( GP) 


i=1 


Where 


Bi» Pis Yi» Ei and Çi are the emission coefficients of generator i 
Piis the power generated by unit i MW 

Nis the number of generating units 

Ezis the emission functieon of unit i 


Eris the total emission, ton/hr. 


The power balance equation[3] is given by 
2, ii Pi=Pp 

Adding loss factor to the equation. 
J$, Pi- Po-Pi 

Power loss PL is calculated as 


PL=) ġa > Pi BiP; + X, Boi Pi+ Boo 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


The actual power generation for the generator will be between its maximum and minimum limits which is 


represented as 


pmi» < Pi < pmax 


(10) 
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3. WHALE OPTIMIZATION ALGORITHM 
In The mathematical model of whale optimization [6] is detailed in this section. 


3.1. Inspiration 

The hunting method of humpback whales, known as bubble-net-feeding, is used as an inspiration to 
create this algorithm. This method highlights the intelligence and corporation of humpback whales. The 
humpback whales hunt in groups (the group size may be as high as a dozen of humpback whales) and the big 
whale which is the leader finds the group of fishes which is to be hunted. The group goes underneath the 
water and the lead whale, who is the bubble blower, produces spiral big bubbles towards the surface of the 
water. These spiral bubbles (shape of ‘9°) interrupt the fishes to swim through the same and they got struck 
up within the spiral bubble. Now the group of whales comes out of the water with the open mouth towards 
the water surface, within the spiral bubble, moving synchronously and hunting all the fishes within the spiral 
bubble. All the whales forms exactly the same position with respect to the lead whale while hunting. The 
bubble-net feeding is a unique behavior that can only be observed in humpback whales. In this work, the 
spiral bubble-net feeding maneuver is mathematically modeled in order to perform optimization. 


3.2. Mathematical model and optimization algorithm 

In this section, the mathematical model of complete behavior of humpback-whale-hunting is done 
which includes search for small fishes, encircling prey and spiral bubble-net feeding maneuver. The Whale 
Optimization Algorithm (WOA) is then proposed. 


3.2.1. Encircling prey 

Humpback whales can recognize the location of prey and encircle them. Since the position of the 
optimal design in the search space is not known a priori, the WOA algorithm assumes that the current best 
candidate solution is the target prey or is close to the optimum. After the best search agent is defined, the 
other search agents will hence try to update their positions towards the best search agent. This behavior is 
represented by the followingequations: 


D = |CX*(t) —X(t)| (11) 


X(t +1) =X*(t) -A-D 
(12) 


where 


‘t’ indicates the current iteration, 

‘A’ and ‘C’ are coefficient vectors, 

X* is the position vector of the best solution obtained so far, 

X is the position vector, | | is the absolute value, and - is an element-by-element multiplication. 
X» should be updated in each iteration if there is a better solution. 


The vectors A and C are calculated as follows: 


A=2d-r-d@ 
(13) 


(14) 


where 4 is linearly decreased from 2 to 0 over the course of iterations (in both exploration and exploitation 


phases) and d is a random vector in [0,1]. 

Figure 1(a) illustrates the rationale behind Equation (12) for a 2D problem. The position (X,Y) of a 
search agent can be updated according to the position of the current best record (X*,Y*). Different places 
around the best agent can be achieved with respect to the current position by adjusting the value of A and C 
vectors. The possible updating position of a search agent in 3D space is also depicted in Figure 1(b). It should 
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be noted that by defining the random vector (7 ) it is possible to reach any position in the search space 
located between the key-points shown in Figure 2. Therefore, Equation (12) allows any search agent to 
update its position about the current best solution and simulates encircling the prey. 

The same concept can be extended to a search space with n dimensions, and the search agents will 
move in hyper-cubes around the best solution obtained so far. As mentioned in the previous section, the 
humpback whales also attack the prey with the bubble-net strategy. This method is mathematically 
formulated in the next section. 


(X*-AX.Y*) e (X.Y*) 
104 


(X*-AX, ¥*-AY) (X*,Y*-AY) (XY*-AY) 


Figure 1. Bubble-net search mechanism implemented in WOA (Xx is the best solution obtained so far): 
(a) shrinking encircling mechanism and (b) spiral updating position 


3.2.2. Bubble-net attacking method (exploitation phase) 
To mathematically model the bubble-net behavior of humpback whales, two approaches are 
designed as follows: 


Shrinking encircling mechanism: This behavior is achieved by decreasing the value of T in the 


Equation (13). Note that the fluctuation range of A is also decreased by © . In other words A is a random 
value in the interval [—a,a] where a is decreased from 2 to 0 over the course of iterations.Setting random 


values for A in [—1,1], the new position of a search agent can be defined anywhere in between the original 
position of the agent and the position of the current best agent. Fig. 3(a) shows the possible positions from 
(X,Y) towards (X*, Y*) that can be achieved by 0< A <1 ina 2D space. 

Spiral updating position: As can be seen in Figure 3(b), this approach first calculates the distance 
between the whale located at (X,Y) and prey located at (X*,Y*). A spiral equation is then created between 
the position of whale and prey to mimic the helix-shaped movement of humpback whales as follows: 


X(t +1) =D -e . cos (271) +X (t) (15) 


Where D' = |X* (t) -X (6)| and indicates the distance of the i th whale to the prey (best solution obtained 
so far), b is a constant for defining the shape of the logarithmic spiral, | is a random number in [—1,1], and ‘.’ 
is an element-by-element multiplication. The humpback whales swim around the prey within a shrinking 
circle and along a spiral-shaped path simultaneously. To model this simultaneous behavior, we assume that 
there is a probability of 50% to choose between either the shrinking encircling mechanism or the spiral model 
to update the position of whales during optimization. The mathematical model is as follows: 


X*(t) -A-D if p<05 


X(t+1)={_, = 
D' -e".cos(2ml)+X*(t) if p>0.5 


(16) 


where p is a random number in [0,1]. In addition to the bubble-net method, the humpback whales search for 
prey randomly. The mathematical model of the search is given in the next section. 
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3.2.3. Search for prey (exploration phase) 

The same approach based on the variation of the A vector can be utilized to search for prey 
(exploration). In fact, humpback whales search randomly according to the position of each other. Therefore, 
we use A with the random values greater than 1 or less than -1 to force search agent to move far away from a 
reference whale. In contrast to the exploitation phase, we update the position of a search agent in the 
exploration phase according to a randomly chosen search agent instead of the best search agent found so far. 
This mechanism and |A| > 1 emphasize exploration and allow the WOA algorithm to perform a global search. 
The mathematical model is as follows: 


(17) 


X(t +1) = Xna — A.D (18) 


4. WAO FOR ECONOMIC DISPATCH CASE 

Merit order dispatch case is implemented using WOA algorithm. The Economic dispatch case 
implementation is done on standard IEEE 30 bus system. As it’s a standard test system, various parameters 
had already been recorded. This system was used in many comparable studies in the merit order dispatch. 


4.1. IEEE-30 bus system - results 

This system has 6 generator buses at bus 1, 2, 5, 8, 11 and 13. Each of these generators has their 
own fuel and emission coefficients. These are represented in per unit values to simplify calculations. The 
base considered is 100 MVA. The total demand considered is 2.38 p.u. The various generation parameters are 
tabulated as follows, Generator cost coefficients for IEEE 30 bus system is provided in Table 1 given below. 


Table 1. Sample Cost Coefficients 
Unit æ b c e F 


Pimin Pimax 


1 10 200 100 15 6.283 0.05 0.5 
2 10 150 120 10 8.976 0.05 0.6 
3 20 180 40 10 14.784 0.05 1 
4 10 100 60 5 20.944 0.05 1.2 
5 20 180 40 5 25.133 0.05 1 
6 10 150 100 5 18.48 0.05 0.6 


Generator emission coefficients for IEEE-30-bus system is provided in Table 2 given below. 


Table 2. Sample Emission Coefficients 


Unit a £ Y E 

1 4.091 -5.554 6.49 2.00E-03 2.857 
2 2.543 -6.047 5.638 5.00E-04 3.333 
3 4.258 -5.094 4.586 1.00E-06 8 

4 5.326 -3.35 3.38 2.00E-03 2 

5 4.258 -5.094 4.586 1.00E-06 8 

6 6.131 -5.555 5.151 1.00E-05 6.667 


TEEE-30 bus system simulation and implementation of Economic dispatch using WOA is done 
using MATLAB program. The following are the results obtained in comparison with PSO algorithm. 

The value is found at each step and results are calculated varying the no of iterations and graph is 
plotted for the same. The best result is obtained at 106" iterations and there is no difference afterwards. 
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Table 3. Economic Dispatch results using WOA 


PSO WOA 
Pi 0.099441 0.098762 
Pa 0.36248 0.37468 
P; 0.48349 0.47011 
P, 0.87359 0.92773 
P 0.66428 0.70221 
P; 0.39004 0.38991 
Total Fuel Cost 

40752.19 40232.26 
EEr) INR/hr 
Total Emission 
E 0.213921 0.213841 
(“T ) ton/hr 
Total Emission Cost 

661.018 660.7687 
ECT) INR/hr 
Total Cost TC; 41413.20 40893.03 
INR/hr 


INFERENCES 
Whale optimization provides excellent results for an economic dispatch problem in terms of cost 


optimization and easy convergence. On comparison, the cost optimization of WOA is better than PSO by 
1.3%. There can be further improvements on the algorithm considering the following factors. Use 
combinational algorithms — Suggested algorithms are ALO and WOA extracting effective features of both. 
Include other characteristics of humpback whales on this algorithm. This includes effects of ‘double loop’ 
bubbles into this algorithm. 
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